Mapping differential reddening in the inner Galactic globular 

cluster system ^ 

Javier Alonso- Garcia 

Departamento de Astronomia y Astrofisica, Pontificia Universidad Catolica de Chile, 

782-0436 Macul, Santiago, Chile 
Department of Astronomy, University of Michigan, Ann Arbor, MI 48109-1090 

j alonsoOastro . puc . cl 

Mario Mateo 

Department of Astronomy, University of Michigan, Ann Arbor, MI 48109-1090 

mmateoQumich . edu 
Bodhisattva Sen 

Department of Statistics, Columbia University, New York, NY 10027 

bodhiOstat . Columbia. edu 

Moulinath Banerjee 
Department of Statistics, University of Michigan, Ann Arbor, MI 48109-1107 

moulibOumich . edu 

Kaspar von Braun 

NASA Exoplanet Science Institute, California Institute of Technology, Pasadena, CA 

91125-2200 

kasparOcaltech . edu 



Received 



accepted 



- 2 - 



ABSTRACT 



A serious limitation in the study of many globular clusters - especially those 
located near the Galactic Center - has been the existence of large and differen- 
tial extinction by foreground dust. In a series of papers we intend to map the 
differential extinction and remove its effects, using a new dereddening technique, 
in a sample of clusters in the direction of the inner Galaxy, observed using the 
Magellan 6.5m telescope and the Hubble Space Telescope. These observations 
and their analysis will let us produce high quality color-magnitude diagrams of 
these poorly studied clusters that will allow us to determine these clusters' rela- 
tive ages, distances and chemistry and to address important questions about the 
formation and the evolution of the inner Galaxy. We also intend to use the maps 
of the differential extinction to sample and characterize the interstellar medium 
along the numerous low latitude lines of sight where the clusters in our sample 
lie. In this first paper we describe in detail our dereddening method along with 
the powerful statistics tools that allow us to apply it, and we show the kind of 
results that we can expect, applying the method to M62, one of the clusters in 
our sample. The width of the main sequence and lower red giant branch narrows 
by a factor of 2 after applying our dereddening technique, which will significantly 
help to constrain the age, distance, and metallicity of the cluster. 

Subject headings: Globular clusters: general - Galaxy: bulge - dust, extinction - 
methods: statistical - Globular Clusters: individual (NGC 6266 - M 62) 



- 3 - 



Introduction 



The age, chemical and kinematic distributions of Galactic stellar populations provide 
powerful constraints on models of the formation and evolution of the Milky Way. The 
Galactic globular clusters (GGC) constitute an especially useful case because the stars 
within individual clusters are coeval and spatially distinct. But uncertainties in the 
determination of distances to GGCs, along with uncertainties in the determination of their 
physical characteristics such as age or metallicity ca n result when the redd enig produced by 



20051 ) ■ In this context, 



interstellar dust is not properly taken into account (iCalamida et al.l 
differential reddening across the field of the cluster has proven to be difficult to map, and 
the properties of GGCs that suffer high and patchy extinction, e specially those in t he 



direction of the inner Galaxy, have not been accurately measured (jValenti et al. 



20071 ). 



Several a uthors have used differe nt methods to try to map differential extinction in 



GGCs before: 



Piersimoni et al 



the reddening maps; 



f|2002[) have used colors of v ariable RRLyrae stars to c reate 



Melbourne fc Guhathakurtal (120041 ) and 



Heitsch fc Richtlei] (Il999h hav e 



used photometric studies of the stars in the horizontal branch (HB) 



von Braun &: Mateo 



(1200 ih and 



Piotto et al. 



(Il999b| ) have used photometric studies of main sequence (MS), 



subgiant branch (SGB), and red giant branch (RGB) stars. 

In this paper we describe a new dereddening technique, based on this last approach, 
but with the improvement of calculating the extinction not on a predetermined grid but 



^Based partly on observations with the NASA/ESA Hubble Space Telescope, obtained at 
the Space Telescope Science Institute, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under NASA contract NAS 5-26555. This paper also includes 
data gathered with the 6.5 meter Magellan Telescopes located at Las Campanas Observatory, 
Chile. 
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on a star by star basis. In section 2 we extensively explain this dereddening method and 
in section 3 we apply it to an example cluster, M62. In the next papers of this series 
we plan to apply this dereddening method to a sample of clusters in the direction of the 
inner Galaxy that in principle can be affected by this effect, and analyze the obtained 
results. In Paper II we will present a new photometric database consisting of a sample of 
25 inner GGCs observed using the Magellan 6.5m telescope and the Hubble Space Telescope 
{HST), describe the processes we follow to obtain an accurate astrometric and photometric 
calibration of the stars in these clusters, and, after applying the dereddening technique, 
give the extinction maps along the field of the clusters in the sample, along with new 
cleaner, differentially dereddened color-magnitude diagrams (CMD) of these GGGs. In 
Paper III we will provide an analysis of the stellar populations of these clusters based on 
their dereddened CMDs. Finally, in Paper IV, we will characterize the interstellar medium 
along the low-latitude lines of sight were the clusters lie, based on the properties of the 
extinction maps we derived. 



2. Dereddening technique 

Our dereddening technique is composed of five iterative steps (see figured]): 

• We assign a probability to the stars in our observed regions to belong to the cluster 
or to the Galactic field depending on their positions in the sky and in the CMD (see 
section 12. ip . 

• We build a ridgeline for the stellar population of the observed cluster. The ridgeline 
represents the evolutionary path that a given star is going to follow in the CMD (see 
section 12. 2p . 

• We assign an individual extinction value to every star, based on its displacement from 
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the ridgeline along the reddening vector (see section [273]) . 

We smooth the different individual color excesses over the observed field to generate 
an extinction map (see section [2^ . 

We apply the extinction values from the map to the stars in our observations to 
construct a dereddened CMD (see section 12.51) . 



We repeat these steps iteratively. Every iteration the ridgeline is more accurately defined, 
and therefore the extinction map and the dereddened CMD are also more precise. We end 
up the process when there is a convergence in the calculated ridgeline. 



This technique is base d on work by 



von Braun fc Mated (|200l|) and 



(ll999bl ). and pioneered by iKaluzny fc Krzeminskil ( 119931 ). but unlike them, we do not 



Piotto et al. 



divide the field in a grid of well established subregions a priori, but use a non-parametric 
approximation to smooth the information about the reddening in the field provided by 
every star, without such hard edges. 

The success of this analysis depends critically on the assumption that the stellar 
populations are uniform within individual GCs. Recent studies tend to suggest that this is 
not strictly true, especially in the most massive GGCs and GGC s with an extreme blue 



horizontal branch (EBHB) (IBedin et al 



2004 : 



Piotto et al. 



20071 ). In the cases discovered. 



it is conjectured that the spread in helium content betwe en stars in the cluster is i mportant 



20081 ). making 



but the spread in age and metallicity seems to be small ( iD'Antona fc Caloil 
the spread of the stars in the SGB and upper RGB regions of the CMD also small. Since 
most of the information in our method comes from stars in these regions of the CMDs, 
our technique should not be significantly affected. What is more, in our method variations 
in the reddening must be spatially related, which should not be the case if there are 
multiple helium enriched populations within the cluster. Because of these reasons and 
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since our sample only includes a few massive and EBHB clusters, we expect the differential 
populations effects to be comparatively minor and not significantly affect our dereddening 
approach. 

Also it is important to reahze that our dereddening method does not estabhsh the 
absolute extinction toward a target cluster, so to estimate the absolute extinction in each 
case we need to use other methods, e.g., comparison of our differentially dereddened CMDs 
with isochrone models, or use of reddening estimates from RR Lyr stars in our fields. 



2.1. Field-cluster probability assignment 

Prom our photometric studies we can assign probabilities to the stars in our observations 
to belong to the clusters, based on r, their position in the sky with respect to the center 
of the cluster, and also on (c, m), their color and magnitude position on the CMD. These 
probabilities will be used in the next steps of our technique to give higher weight to the 
information provided by stars with high probabilities of being cluster members. 

Before going any further, it is convenient to explain the notation that we are going 
to use, and to remember the basic rules of probability calculation. The marginal or 
unconditional probability of an event A happening is expressed P{A). The joint probability 
of an event A and an event B happening at the same time is expressed P{A, B). Finally, 
the conditional probability of an event A happening given the occurrence of some other 
event B is expressed P{A\B), and can be written as 

m|B) = :M (1) 

or, using Bayes theorem, as 

. Emm (,) 

Also, notice that in many parts of our analysis it is going to be more convenient to use 
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probability density functions, p{x), to express the probability. These functions describe the 
relative likelihood for a random variable X to occur at a given point in the sampled space 
Q and satisfy 

P{a<X<b)= I p{x)dx (3) 

J a 

For all X : p(x) > (4) 
p{x)dx = 1 (5) 
Whenever needed in this work, th e probability density functions are calculated non- 



parametrically using locfit. Locfit (ILoader 



19991 ) is a local likelihood estimation software 



implemented in the R statistical programming language (see appendix). Locfit does not 
constrain the probability density functions globally, i.e., it is non-parametric, but assumes 
that locally, in a local window around a certain point in the sample space, the functions can 
be well approximated by a polynomial. Hence the main inputs that locfit requires are the 
positions of the observed elements in the sample space, and a parameter to define the local 
window where the function is approximated parametrically, i.e., a smoothing factor given 
by the maximum of two elements: a bandwidth generated by a nearest neighbor fraction, 
and a constant bandwidth. 

Now, if we suppose the variable X indicates membership to the cluster {X = 1 
indicates observation of a member, and X = indicates observation of a non- member), in 
the next paragraphs we calculate P{X = l\r,c,m), the conditional probability of the stars 
being members of the cluster, given their position in the sky and in the CMD. Actually, it 
is more convenient for us to calculate first P{X = l|r) and P{X = l\c,m), the conditional 
probabilities of the stars to belong to the cluster given just their position in the sky, and 
just their position in the CMD, and only after attempt to calculate P{X = l\r,c,m). 

The conditional probability of the stars being members of the cluster, given just their 
position in the sky as a function of the distance to the center of the cluster, P{X = l|r). 
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can be written, using Bayes theorem, as 

p{r) 

For any given cluster, both p{r\X = 1) and P{X = 1) are unknown. But since P{X = 1) is 
the ratio of the observed cluster stars to the total observed stars, we can rewrite eg. 6 as 

P(X = llr) = ^^^^^ (7) 
St{r) 

where Smem{f) is the number of observed member stars as a function of the distance to the 
cluster center, and St{r) is the total number of observed stars as a function of the distance 
to the cluster center. We still do not know Smem{r), but if we divide both numerator and 
denominator in eg. 7 by A{r), the area coverage at a given r of the field of view (FOV) in 
our observation of the cluster, we can rewrite eg. 7 as 

P(X = l|r) = ^77l''^ (8) 
ft{r) 

where fmem{f ) is the surface density of observed member stars as a function of the distance 
to the cluster center, and ft{r) is the total surface density of observed stars as a function o f 



19621) 



the distance to the cluster center. Now, if we assume an empirical King profile (iKing 
for the cluster, we can write the surface density of member stars as 

I if r > Tf 

where Vc and rt are the core radius and the tidal radius of the cluster. If we make another 
assumption and consider a constant surface density of stars for the non-member population, 
which is reasonable due to the small size of our field, then 

fnonir) = C (10) 

and we have the following functional form for the total surface density distribution of stars 
observed in our FOV 

ftir) = kK{r) + c (11) 
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and, according to eg. 8, we can write the conditional probability of the stars being members 
of the cluster, given their position in the sky, as 



P(X = = (12) 

kK{rj + c 



or, alternatively. 



P(X = l|r) = 1 - P(X = 0|r) = 1 - (13) 

kK[r) + c 

Therefore to calculate P{X = l|r) we need to find the coefficients k and c, since K{r) can 
be obtained ta king the valu es for rc and rt for every cluster in our sample provided in the 



Harris catalog ( iHarris 



19961 ) ■ These two remaining coefficients, k and c, can be easily found 
by a least squares fit of K{r) to the observed surface density of stars ft{r) (see eg. 11, 
figure |2l and tabled]). To find the observed ft{r), we need to first calculate the number of 
stars St{r), and then divide it by the area coverage A{r). As we mentioned when passing 
from eg. 6 to eg. 7, the number of stars St{r) is just the probability density p(r) of stars 
in our observation, multiplied by the total number of observed stars. We calculated the 
probability density p(r) of stars by feeding locfit with the radial positions of the stars and a 
smoothing factor with a constant bandwidth of 0.25 arcmin. The surface coverage A{r) in 
our observation is not trivial to find, since our FOV is a square not centered in the cluster 
center. In order to quickly calculate this area, we create a grid of points equally spaced over 
our FOV, calculate the probability density p(r) of points, and then multiply by the area 
of the FOV. We calculated the probability density p(r) of points by feeding locfit with the 
radial positions of the points in the grid and a smoothing factor with a constant bandwidth 
of 0.25 arcmin. 

The conditional probability of the stars being members of the cluster, given just their 
position in the color magnitude diagram, P{X = l|c, m), can be written, using Bayes 
theorem, as 

p(c,m) 
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For any given cluster, we do not know the true distribution in magnitude and color in 
our observations of genuine GC members, but we can model the distribution of the field, 
non-member stars. We use the model of the Galaxy described in 



Robin et al. 



(120031), from 



now on referred as the 'Besangon model', to obtain p{c,m\X = 0). We can now rewrite 
eg. 14 as 

p{c,m\X = 0)P(X = 0) 



P{X = l|c, m) = 1 - P{X = 0|c, m) = 1 



p{c,m) 



(15) 



where p{c, m|X = 0) is obtained feeding locfit with the position in the CMD of the Galactic 
field stars provided by the Besangon modej^ and a smoothing factor with a nearest neighbor 
ratio of 0.01, p(c, m) is obtained feeding locfit with the position in the CMD of the stars in 
our observations and a smoothing factor with a nearest neighbor ratio of 0.01, and P{X = 0) 
is just the ratio of the total number of modeled non-member stars for an area equal to the 
FOV of the observation, to the total number of observed stars, which include field and 
cluster stars. To take care of the different variation ranges in color and in magnitude for 
the stars in the CMD a scale (see appendix) of 1 in color to 5 in magnitude is also provided 
to locfit as an input parameter to calculate p{c,m\X = 0) and p{c,m). In figure [3] we can 



Similar methods to calcu 
ba sed on calculation s by 



by 



Law et al. 



^1^, ..^J ^^^^ ^^^^^ ^^^v. v^^ v-^-v^ -^^^^ 

ate P(X = 1 c, m) have been employed by ^ 


Huffhes et al. (2007) 


Hughes & Wallerstei 


n f2000) and Miffhell et al. (1998 


), and also 


d on calculations by 


Odenkirchen et al. (200! 


.) and 


Grillmair et al. 



( 119951 ). In contrast to these other studies we do not eliminate stars from our analysis based 
on these probabilities, but downweight them. Also, as we show at the end of this section, 
we do not just use P{X = l|c, m), but P{X = l|r, c, m) in our following analysis. 



Given the importance of the Besangon model (and the CMDs it provides) in analyzing 



^These CMDs can be easily obtained via the web interface provided at: 
http:/ /model.obs-besancon.fr/ 
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the Galactic field component of our observations and calculating p(c, m\X = 0), we describe 
here in more detail the input parameters that we use in the web interface of the Besangon 
group to model the non-member stellar populations in the vicinity of the sampled clusters. 

• Field of view. We use the small field option from the web interface, in which the 
modeled stars are all supposed to be at the same coordinates (see table |2]), implicitly 
assuming that the field star density gradient across the FOV is negligible. The same 
assumption was made in the calculation of fnon{f) (see eg. 10). To obtain higher 
statistical significance in the models, we can increase the selection area to obtain a 
larger number of stars in the model. p(c, m|X = 0) is independent of the area and 
the number of stars used, since it is normalized (see eg. 5). In practice, we set this 
parameter to select ~ 300000 stars (see table [2]). We have to notice though, that 
when calculating P{X = 0) in eg. 15 we need to normalize the number of stars in the 
model to the FOV. 

• Extinction law. We apply locfit to map p(c, m) as a function of c and m using the 
CMDs from our observations. From the modeled stars, we can also use locfit to map 
p(c, m|X = 0) as a function of c and m, realizing that this map changes depending on 
what extinction Ay we assumed for our modeled data. So if we define gext as 



we need to realize that this function depends also on the extinction assumed for the 
modeled CMDs, i.e., gexti^v)- For the fields studied, we choose the modeled stars to 
simulate no interstellar extinction in the beginning. Then we move the density maps 
of the models along the reddening vector in increments of extinction AAy = 0.01, 
calculating g^xt for every increment. The value of Ay where ge^t is a maximum is the 
value we choose for the absolute extinction of our observation. This Ay is generally 




(16) 



c.m 
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equal or within a few hundredths of a magnitude to the one provided in iHarrisI ( 119961 ) . 
The extinction is then simulated in the model as a cloud with an extinction equal to 
this Ay at distance of pc (see table [2]). 

• Characteristics of the stars used in the model. We allow all ages, Galactic 
components, and spectral types provided by the Besangon model. We provide a 
limit only in the intervals of magnitudes, equal to the ones of our observations (table 
[2]). We also model our observational errors Cph with an exponential function of the 
apparent magnitude m in each band: 

apf^ = A + exp{Cm - B) (17) 

so the models will take the photometric error into account providing a more realistic 
approximation (see table [2]). 

Certain problems arise in this procedure of calculating the conditional probabilities 
P{X = l\r) and P{X = l|c, m), including: 

• We have to take into account the completeness factor of our photometric observations. 
In general, completeness decreases with magnitude somewhat smoothly at the faint 
end, but suddenly at the saturation limit. Because of this, and since we are trying to 
calculate the probabilities comparing our observations with models (King profile and 
Besangon model) we need to find where the completeness factor starts to decrease 
at the faint end. Usually this is done with artificial-star tests, injecting stars in the 



CMP and analyzing the amount and magnitudes of the ones recovered (jPiotto et al. 



1999al ). But our large sample and large cluster star density gradients within individual 
fields make artificial star tests highly inefficient, so we explore another method that 
is well-suited to our specific fields, and also much more efficient to implement. From 
the Besangon model we can calculate the number of non-member stars per area in 
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the field that we have in a given magnitude range, and from our data we can also 
find the number of non-member stars per area in the field in a given magnitude 
range just by fitting eg. 11 in that magnitude range. So we iteratively calculate 
both numbers making the magnitude range smaller taking away stars in the 0.1 mag 
fainter end every time we do the calculation, repeating the process until we are able 
to get the cumulative luminosity functions (LF) for the seven fainter magnitudes for 
the non-member stars from every one of our pointings. Then we can compare the 
modeled and observed LFs to see where the functions deviate one from the other. 
But this approximation proves to be not very accurate, especially in cases in which 
the number of non-member stars in the field is small. Instead we explore yet another 
approximation. We notice that the LF in the Besangon model is always a concave 
function, while the observed LF of the non-member stars in the field becomes a convex 
function at the faint end. We can then identify the limit where the completeness 
factor starts to decrease as the inflection point where the observational cumulative LF 
changes from concave to convex. The inflection point is the point where the derivative 
of the LF has a extreme, a minimum in our cases. And the derivative is easily found 
using locfit to calculate non-parametrically the LF and its derivative (see appendix). 
The completeness limit is usually located in our observations at ~ 2 magnitudes from 
the faint end in V (see table [3]). In the next steps of our dereddening method we only 
use stars brighter than this limit. 

The completeness factor also depends on the distance to the center of the cluster, 
strongly so for the more crowded cases. Close to the cluster center we are not only 
missing stars at the faint end of the magnitude distribution, but at all magnitudes, 
since the lower completeness fa ctor in this region is mostly caused by masking 



saturated regions by DoPHOT (ISchechter et al. 



19931 ). the program used to do the 



photometric analysis of the images. From the calculation of the P{X = l|r), we can 
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see that this reduced completeness factor reveals itself as a deviation from the model 
fit shown in eg. 11 at distances close to the center of the cluster (see figure [2]). If we 
assume that we are missing an equal percentage of stars from the cluster and from 
the field, we can calculate its effect as the ratio between the observed and modeled 
ft{r). In the CMD, this lower completeness factor can lead to a miscalculation of 
P{X = l\c,m) if the FOV is small. However, since the area of the regions whose 
extinction we are able to map are much bigger than the areas where the reduced 
completeness factor is an issue, the comparison of our observations with the Besangon 
model is not significantly affected. The only appreciable consequence is that the 
calculated P{X = 0) is a little higher than the actual P{X = 0). But this high value 
for P{X = 0) can be corrected if when we calculate the number of stars, we take 
into account the deviation of the observation from the model fit as a function of the 
distance to the center of the cluster, as we mentioned before. 

P{X = l|r) is obtained usin g the values for the core and tidal radii provided in the 



Harris catalog (IHarris 



19961 ). for every cluster in our sample. Notice that these values 
are not all of equal accuracy. Small variations in the initially adopted and could 
produce small changes in the individual P{X = l\r) values, but they will not alter 
the general trend (stars closer to the center of the cluster have higher P{x = l|r) than 
stars located further away), and therefore they will not significantly affect the result 
of our dereddening technique. 

P{X = l\c,m) is obtained comparing a real observation with a model, where 
photometric errors are incorporated by mimicking the real photometric errors by the 
exponential function described above. In general, we observe that stars from the 
models tend to be a little more concentrated along the different field evolutionary 
sequences in the CMDs than stars from real observations. Therefore, for a given 
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magnitude, p{c,m\X = 0) peaks higher and decreases faster in the model than in the 
real observation, producing two main errors in our calculations of P{X = l|c, m). 
The first error is that we mistakenly assign a probability of being in the field which 
is too low, and therefore a probability of being cluster members which is too high, to 
stars in the regions where the observed color is an extreme for a given magnitude (see 
figure [3]). But since there are not many stars in these regions and they usually have 
higher photometric errors in color than our functions describe, they are downweighted 
or removed from our calculations later on (see section [2^ . More important is the 
second error, the case were p(c, m\X = 0) peaks higher in the model than in reality, 
making the probability of the star to be a non-member in the field too high, and 
P{X = l\c,m) too low, lower than 0. In order to avoid these non-physical cases, and 
even too extreme cases that can give too low a weight value to the information coming 
from those stars, we put the hard limit P{X = l|c, m) > 0.1 (see figure [3]). 

Now we are in conditions to calculate P{X = l|r, c, m), the conditional probability of 
the stars being members of the cluster given both their position in the sky and in the CMD. 
According to Bayes theorem, 

p(r,c,m\X = 0)P(X = 0) , , 

P{X = l|r, c, m) = 1 - P{X = 0|r, c, m) = 1 - ^-^^ — ^ '-^ '- (18) 

p{r,c,m) 

and since the positions in the sky and in the CMD for non-member stars are independent □ 

P(A- = i|..c,™) = i-^<I^^^M£i!!4£^W^ (19) 

p(r, c, m) 

and using Bayes theorem again (see eg. 6) we can write 

p(x = i|.,c,™) = i-^<Hii^^^^Mrm^ (20) 

p(r, c, m) 



^Notice that in the denominator p{r,c,m), the variables are not independent, i.e., the 
position in the CMD for a star in our observation is related to its spatial position in the sky, 
due to the mixture of cluster and field populations. Hence p(r, c, m) 7^ p{r)p{c,m). 
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P{X = 0|r), p(r), and p{c,m\X = 0), the three elements in the numerator have been 
calculated previously, so now we just need to find the denominator p(r, c, m), the probability 
density of stars in our observation as a function of both position in the sky and in the 
CMD. To calculate it we feed locfit with the position in the sky and in the CMD of the 
stars in our observations and a smoothing factor with a nearest neighbor ratio of 0.1. Since 
the units of the three physical magnitudes r,c,m are different, we need to feed locfit with 
a scale factor (see appendix) to calculate the distance between neighbors in the r, c, m 
space. Experimentally we found that a scale of 1 magnitude in color to 5 magnitudes in 
brightness and 10 arcmin in distance to the the GC center provides good results, and it is 
the one we used as an input parameter of locfit to calculate p{r,c,m). The results of the 
P{X = l\r,c,m) calculation are shown in figure HI We can observe that the same problems 
that were previously mentioned affecting the calculation of P{X = l|r) and P{X = l\c,m) 
are present here. But for most of the stars observed, this formalism allows a more precise 
way to discriminate between field and cluster stars than looking to just P{X = l|r) or just 
P{X = l\c,m). 

2.2. Building the ridgeline 

Stars that are still alive in globular clusters evolve along the CMD following a well 
determined path. They spend most of their lives in the MS burning hydrogen in their 
core. They leave the MS at the turn-off point (TO) growing in size and burning hydrogen 
in a shell as they move across the SGB and the RGB. Then they move to the HB where 
they burn helium in the core and hydrogen in a shell. They leave the HB moving across 
the asymptotic giant branch (AGB) burning helium and hydrogen in a shell, then moving 
rapidly across the post AGB phases to finish their lives as white dwarfs (WD). The path of 
this standard evolutionary sequence in the CMD of a globular cluster is shown in figure [51 
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We try to model the first stages (MS, SGB and RGB) of tlie evolutionary patli fol 



owed 
to 



by tlie stars in a GC building the ridgeline for the CMDs of the GGCs, using locfit 
construct a univariate non-parametric regression of the color of the stars as a function of 
the magnitude. This process is composed of three iterations. 

In the first iteration, we carry out a non-parametric regression in the whole CMD to 
obtain a first estimation of the cluster ridgeline as a function of magnitude. We feed locfit 
with the magnitudes and colors of the stars and a smoothing parameter with a constant 
bandwidth of 0.2 magnitudes. Also, we give locfit preliminary weights for the stars 

W = P{X = l\r,c,m)wph (21) 

where P{X = l|r, c, m) is the membership probability calculated in the previous subsection 
and Wph is a photometric weight equal to the inverse of the square of the Poisson error of the 
photometric magnitude Wph = l/cTy. These preliminary weights W are calculated to favor 
the magnitude and color information provided by stars with more accurate photometry and 
higher probability of being cluster members according to both their position in the sky 
and in the CMD. The smoothing parameter values used have proven experimentally to be 
adequate for describing the MS and part of the RGB, but the regression fails to describe 
the SGB, and the stars in the HB create problems when trying to model the whole RGB. 
Still, they produce a reasonable first approximation (see upper left plot in figure [5]) and let 
us identify the TO point, which we define as the bluest point of our ridgeline. We should 
note here that the MS is defined not to the faintest magnitude limit of the data, but to the 
completeness limit defined in section 12.11 

In the second iteration we try to get a less noisy whole ridgeline, along with a better 



'^Locfit also calculates regressions non-parametrically when it is fed with the positions of 
two sets of related observations (see appendix). 
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defined RGB ridgeline. As a first step we divide the initial ridgeline in three regions: the 
MS region, where the ridgehne magnitudes are m > rriTo, the SGB region, where the 
ridgehne magnitudes are itlto — 1 < m < rrixo, and the RGB region, where the ridgehne 
magnitudes are m < ttito — 1- To smooth out the approximation for the ridgeline, we take 
the points from the ridgeline from the first iteration for every region and smooth them 
using locfit with a nearest neighbor ratio of 0.7 {locfit default), no preliminary weights or 
constant bandwidth requirements. The RGB ridgeline is still not well defined, because of 
the effects of stars in the HB. To calculate a better ridgeline for the RGB getting rid of the 
effects of the blue part of the HB, we take only stars in the magnitude range m < rriTo ~ 1 
and in the color range coIto < col. We feed locfit with their magnitudes and colors, and 
a constant bandwidth of 0.2 mag and a nearest neighbor ratio of 0.1 to calculate the 
smoothing parameter. We use the nearest neighbor ratio because if not, the upper part of 
the RGB, usually scarcely populated, can become very noisy. Still, the red part of the HB, 
whenever present, is going to deviate the calculated RGB ridgeline. To eliminate its effect 
we have to locate first where the HB is. We look for a change in the sign of the slope of 
the RGB. The bluest point brighter than the point where the slope sign change is, shows 
where the HB is. And the difference between these two points provides information about 
the thickness of the HB, rriHBthickness- Once the HB is located, we repeat the analysis done 
before for the RGB, with the same smoothing parameter and weights, but now omitting 
points in the interval niHB - rriHBthickness <m< uihe + rriHBthickness- That way we have 
a better ridgeline for the upper and lower regions of the RGB. To get a ridgeline for the 
whole RGB, smoothing any noisy part, we follow a similar process to what we do in the 
beginning of this second iteration. We take the points in the description of the upper and 
lower RGB ridgeline, and use locfit again, but now with a nearest neighbor ratio of 0.7, no 
preliminary weights or constant bandwidth requirements. The ridgeline is much better now, 
but still we are a little off in the description of the SGB (see upper right plot in figure [5]). 
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In the fin al iteration, we try t o calcu late a better ridgeline for the SGB. In order to do 
that we follow 



Marin-Franch et al. 



( 120091 ) ■ where they use rotated histograms to recalculate 
the ridgeline for the GC in their sample, using only stars perpendicular to the ridgeline 
in the calculation. Although our approximation is similar, the implementation is a little 
different. We calculate the slope of the curve at every point of the previously calculated 
ridgeline using locfit to get the derivative of the slope at a given point of the ridgeline. 
Once we have this, we rotate the coordinate system of every star at every magnitude with 
the angle a, of the slope of the ridgeline at that magnitude, and centered at the ridgeline. 
Once this is done, we can get a ridgeline in the new coordinate system using locfit again. 
This new ridgeline should be a straight line along the Y coordinate with a value of 
in the X coordinate. Any deviation from that means that the ridgeline in the original 
color-magnitude coordinate system requires a more accurate calculation. Notice that in 
order for the new coordinate system to consistently show these deviations, we need the 
range of both coordinates to be similar in the old coordinate system. This is not the case 
for the range of colors and magnitudes presented in a CMD. To try to get them to a similar 
scale we multiply the color by a factor of 5 before doing the rotation. In the calculation of 
the ridgeline in the new coordinate system the preliminary weights and smoothing factor 
used for locfit are the same as in the first iteration. We observe that the X coordinate of 
the ridgeline in the new reference system does not deviate significantly from in the region 
of the MS stars, but does so in the SGB and upper RGB region. In the RGB region it is 
expected due to the scarcity of the population there. We took care of that in the previous 
iteration, so we concentrate our attention in the SGB region. We derotate the ridgeline in 
that region (the part where the original magnitudes are in the range rrirpo — 1 < m < rriTo) 
to get their coordinates in the color-magnitude coordinate system, and after smoothing out 
any noise using locfit with a nearest neighbor ratio of 0.7, we put together the different 
parts of the ridgeline (see lower left plot in figure [5]) . The new ridgeline seems to follow 
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accurately the different evolutionary sequences present in the GCs, and serves as the basis 
of our subsequent analysis. 

In section 12.41 we discuss how we can use stars in the HB to test the accuracy of our 
method. To carry this test, we need to model a ridgeline for this region too. The process 
to find the ridgeline here is a little more interactive than for the other regions of the 
CMD. First, we have to decide by visual inspection of the CMD if we are dealing with 
a cluster that has only blue, only red or both sections of the HB. If the cluster shows 
only a blue HB, we carry out a non-parametric regression on the stars bluer than the TO 
to obtain an estimation of the cluster HB ridgeline as a function of magnitude. We feed 
locfit with the magnitudes and colors of the HB, and a constant bandwidth of 0.2 mag 
and a nearest neighbor ratio of 0.1 to calculate the smoothing parameter. If only a red 
HB is present instead, the non-parametric regression is performed on stars in the interval 
mHB - rriHBthickness < m < niHB + mHBthickness to obtaiu an estimation of the cluster 
HB ridgeline as a function of color. The smoothing parameter is calculated by locfit with 
the same parameters as for the blue HB. Finally, if the HB shows red and blue sections, 
we calculate both independently following the methods previously described, and then we 
smooth the result as a function of color, feeding locfit with the magnitudes and colors of 
the points in the description of the blue and red HB ridgeline, and a nearest neighbor ratio 
of 0.25, no preliminary weights or constant bandwidth requirements (see lower right plot in 
figure [5]) . 

2.3. Calculating an extinction for every star 

To calculate an extinction for every star we move the stars along the reddening vector 
until they intersect the ridgehne at an astrophysically reasonable location (see figure [6]). 
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The reddening vectors are described by the equations 



Av/E{B-V) = 3.317 



(22) 



and 



Av/E{V -I) = 2.411 



(23) 



as given in iSch 



CardeUi et ah 



egel et al.i 



If 



(11998 



, wh i ch are evaluated using the Ry = 3.1 extinction laws of 



fll989h and lO'Donnelll ( 119941 ). 



The HB ridgeline, and therefore also stars in the range — ^HBthickness < m < 
TnuB + nT'HBthickness, are not used in this calculation, since we want to use them as an 
independent test of the accuracy of our method (see section [23]). We do not use stars 
dimmer than the completeness limit defined in section [TT] either. These stars and stars that 
do not intersect with the ridgeline are given a weight in the calculation of the extinction 
map described in next section. 



fl200lh 



b calculate the error of these shifts we follow the analysis in ivon Braun fc Mateo 



We created an error ellipse for every star defined by the Poisson error in its color, cTc, 
and magnitude, am- 



Since the color and magnitude errors are correlated, the error ellipse is tilted. The tilt 
angle and the length of the semi-major and semi-minor axes of the error ellipse are 
functions of the error in color and in magnitude. 

2(T^rr7, 



tan(2^) 



+ ^m 



+ 



nl/2 



^ cm 



4 



1/2 



(24) 

(25) 
(26) 
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where 6 is the tilt angle of the error ellipse, Uc' and am' are the semi-major and 
semi-minor axis of the tilted ellipse, and a cm is the covariance of the color and 
magnitude. 

• For every star, we move the error ellipse along the reddening vector. The point of first 
contact (pfc), i.e., the point where the ellipse touches the ridgeline for the first time, 
and the point of last contact (pic), where the ellipse touches the ridgeline for the last 
time, represent the one-sigma-deviation points (see figure E]). 

• Although these two contact points are not necessarily symmetric about the reddening 
value of the star, i.e., the center of the error ellipse, we are going to define the error 
in the shift as: 

(^ellipse = 0.5(p/c + plc) - center (27) 

2.4. Creating the extinction map 

We now take all the color excesses for the individual stars and smooth them using 
locfit to build a bivariate non-parametric regression of the extinction as a function of the 
spatial coordinates right ascension and declination, up to a distance from the center of the 
cluster equal to where P{X = l|r) = 0.1 or the limit of our observations, whichever comes 
first (see table |3]). This limit is adopted after preliminary application of our method to the 
clusters in our sample with looser restrictions, and not observing an improvement in the 
dereddened CMDs (see section 12. 5p for stars located beyond these ranges. Stars that were 
given zero weight in the previous subsections (stars in the HB region, stars dimmer than the 
completeness limit, stars that do not intersect with the ridgeline after being moved along 
the reddening vector) are not taken into account for the calculation of the map, although 
after the map is built a reddening correction is applied to all of them. Preliminary weights 
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are assigned in locfit to all stars as they were calculated in previous subsections, but now 
Wph = ^/o'eiiipse- The spatial smoothing kernel in this case has a constant bandwidth of 
0.5' X 0.5' or a nearest neighbor ratio of 0.03, whichever is bigger. These values for the 
smoothing parameter are used after experimentally checking in a few of our clusters which 
kernel provides us with a tighter HB, in the sense of a smaller standard deviation of the 
stars from the ridgeline of this region (see section 12121) . after generating the dereddened 
CMD (see section [231) . This is an independent test to assess the quality of our dereddening 
method, since our technique does not use stars in the HB to calculate the extinction map. 
The tests gave similar results for a range of parameters. Since the tests were not highly 
conclusive on choosing a particular set of parameters, the set used is an average of the best 
values obtained for the different clusters tested. 

Along with the extinction map, an extinction precision map is also calculated. This 
is done assuming that the precision of the extinction map is a function of the spatial 
coordinates. Hence, we can write 

Yi = r(«,. Si) + aiai, 5,)ei, E(q) = 0, V(q) = 1 (28) 

where Yi is the calculated extinction for the i-th star and r{ai, Si) is the value for the 
extinction provided by the regression function at the coordinates {a, 6) of the i-th star. 
Therefore, regressing directly the square residuals {Yi — r(aj,5j))^ provides us with an 
estimate a'^{ai,Si) of the variance, and its square root allows us to plot a map of the 
standard deviation a as a function of the spatial coordinates. This precision map provides 
us with a tool to improve our extinction map. We iteratively repeat the whole process of 
the calculation of extinction maps, but every time using only stars that have extinctions 
that are no more than 3a away from the value that our previous extinction map gives for 
those coordinates. Preliminary weights and smoothing parameters are the same in every 
calculation. We repeat the process until it converges, i.e., until the number of stars used in 
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the calculation of the extinction maps changes by less than 1% with respect to the previous 
iteration. This way for every cluster we are able to provide an extinction map (see top plot 
of figure [7]), a precision map for this extinction map (see bottom plot of figure [7]), and a 
resolution map with the bandwidth that we have use in a certain region to calculate the 
extinction map there (see central plot of figure [7]) . 



2.5. Creating the dereddened CMD 

All the observed stars are given a reddening correction from the extinction map and a 
dereddened CMD is constructed (see figure [8]). This dereddened photometry is the input 
on the next iteration to calculate the cluster membership probabilities, and to build a new 
improved ridgeline (see figured]). Notice though that the input on the calculation of the 
individual reddenings to generate the extinction maps is the original photometry in every 
iteration (see figured]). 

We iteratively repeat this process until there is a convergence in the calculated ridgeline 
(see figured]), which results in good convergence of the reddening values of the extinction 
map. For our dataset, convergence usually occurs after just two or three iterations. 



3. Example: Dereddening M62 

The method will be extensively applied to a sample of 25 globular clusters in following 
papers. Here, as an example, we applied the method to one of those clusters, M62. 



M 62 (NGC 6266) is a mod erately redde ned, E{B - V) = 0.47 clus ter, located only 



1.7kpc from the Galactic center (iHarris 



2005 



19961 ). with a patchy extinction (Contreras et al. 



2010h . The observation, reduction and calibration processes will be explain in detail 
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in Paper II. Suffice it to say here that B, V, and / data for this cluster where obtained 
with the IMACS camera in the Magellan 6.5m telescope, and, in V and I for the inner 
part of the cluster with the WFPC2 camera on the HST. After applying fairly standard 
reduction and calibration processes to these data, we merged the databases from the two 
telescopes to obtain accurate astrometry (cr ~ 0.25" dispersion) and accurate photometry 
(a ~ 0.2 magnitude dispersion) in B,V and I for the stars in a field of 15' x 15' centered in 
the cluster. In figure [8] we can see the resulting CMDs of this cluster. 

The ffist step of our dereddening method is to calculate the probabilities P{X = 
l\r,c, m) for the observed stars to belong to the cluster. To do that, we use the probability 
densities of the stars as functions of distance to the center of the cluster and of color and 
magnitude, a King proffie plus constant field model, and a Besangon model for the field, as 
explained in section \27L\ (see figures El [3l andHj). The different parameters used to build 
the King and Besangon models for this calculation are shown in tables [1] and [2l Notice 
that we restrict our analysis to stars in the limits shown in table [3l because of the reasons 
mentioned in the previous sections. 

Once the probabilities to belong to the cluster have been calculated, we need to build 
the ridgeline. Following the the three step recipe explained in section 12. 2[ we get the 
ridgelines in the two available colors (see figure [5]). After that, we move the stars along the 
reddening vector until they intersect the ridgeline (see figure [6]) and smooth the resulting 
individual color excesses, as explained in sections 12.31 and 12. 4[ We obtain the extinction 
map shown in figure [71 From this map, we take the relative extinction that corresponds to 
every star observed. We then plot the CMD of the stars in our observation after having 
been corrected for the differential reddening and use it as the input for the next iteration. 
The process converges for both colors after 3 iterations and we generate the dereddened 
CMDs that we can see in figure O 
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We leave the detailed examination and analysis of the physical characteristics of the 
cluster to Paper III, but we mention here how much better the definition of the different 
parts of the cluster is. The width of the MS narrows by a factor of 2 after being dereddened 
(see figure [9]). This improvement will help us to obtain better constraints on the cluster age 
and distance when compared to theoretical isochrone models. The RGB is also narrowed by 
a factor of 2 (see figure [H]), which will help when trying to better constrain the metallicities. 

As a reference, we compare the final reddening map obtained using our technique with 
maps and extinction values within the covered region available in the literature (see figure 
[TUl) . We found a general agreement in the identification of regions of h igh and low extinctio n, 



Schlegel et al. 



(119981), 



although our map shows sharper features than the one presented by 
perhaps the result of their mapping the interstellar extinction in the foreground and also in 
the background of the cluster (upper plots in figure [ID]), and smoother values than those by 



Contreras et al. 



(|2010[ ). result of our spatial smoothing of the extinction values (lower plots 
in figure [To]). As we have already mentioned, the technique explained in the paper does not 
calculate the absolute extinction toward a target cluster. Comparisons like the ones shown 
in figure [10] will allow to find the extinction of the adopted reddening zero point in our 
dereddening method. 
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provided by MIDEPLAN's Programa Inicativa Cientffica Milenio through grant P07-021-F, 
awarded to The Milky Way Millennium Nucleus. 
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A. Non-parametric analysis using locfit 

In our work we have extensively used non-parametric statistics to estimate probability 
density and regression functions directly, without reference to a specific form. This differs 
from other approaches in which probability density and regression functions are expressed 
in a parametric way, where the function used to describe them can be written as a 
mathematical formula which is fully described by a finite set of parameters that we have to 
find. 

In our analysis we have used locfit, a local likelihood estimation softwar e implenaented 



in the R statistical programming language. Locfit is extensively explained in iLoaderl (Il999l ). 
so here we just want to give a general idea of how it works and the main tuning parameters 
that we have chosen in our analysis. Locfit does not constrain the functions globally, i.e., it 
is non-parametric, but assumes that locally, around a certain point x, the function can be 
well approximated by a member of a simple class of parametric functions. Locfit defines a 
local window around a point, weighting the observations according to their distance to that 
point 

, , I if x,<\x + hix)\ 

Wi[X) = < 

I if Xj > |x -|- h{x)\ 

and inside this local window, the function is approximated by a polynomial, using not the 
usual local least square criterion but a local log likelihood criterion. We have to choose the 
weighting function and the order of the polynomial. For both cases we take the default given 
by the program: a polynomial of order 2 with a tricube weight function W{u) = (1 — ImI^)^. 
Still we are left with one last argument to choose, the smoothing parameter that controls the 
bandwidth h{x). This smoothing parameter is defined by the maximum of two elements: a 
bandwidth generated by a nearest neighbor fraction < a < 1, and a constant bandwidth. 
The nearest neighbor bandwidth is computed in two steps, first computing the distances 
\x — Xi\ for all the data points and then choosing h{x) to be the hh smallest 
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distance, where k = [na]. 

Also the hkehhood criterion can be chosen, and we choose for the regression the family 
qrgauss, which is equivalent to a local robust least squares criterion where outliers ar e 



iteratively identified and downweighted, similar to the lowess method (I Cleveland 



19791). 



If we want to give more weight to some points than others, locfit allows preliminary 
weights to be given to all observations. 

A scale factor can be applied to the different variables in multivariate fitting, when 
variables are measured in non-comparable units, or when we want to give more importance 
to one of them like in the determination of the densities of stars in our CMD, where the 
range in color for the stars is smaller than the range in magnitudes. 

In addition to calculating the non-parametric function that fits the data, locfit can also 
calculate the derivative of that function. 

Finally, a word is in order about the evaluation structures in locfit. Locfit does not 
perform local regression directly in every point, but selects a set of evaluation points 
obtaining the fit there and interpolating later elsewhere. This is done for efficiency: it is 
much faster to evaluate the structure in a small number of points and then interpolate 
for the rest. But we should make sure that we do not lose information on the process. In 
order to achieve this, locfit uses, by default, a growing adaptive tree, which is the evaluation 
structure that we use in our analysis. A growing adaptive tree is a grid of points. One 
begins by bounding the data in a rectangular box and evaluating the fit at the vertices of 
the box. One then recursively splits the box into two pieces, then each subbox into two 
pieces and so on. For this kind of structure, an edge is always split at the midpoint and 
the decision to split an edge is based solely on the bandwidths at the two ends of the edge, 
depending on the score 5ij = dij/min{hi,hj) where dij is the distance between the two 
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vertices of an edge, and hi and hj are the bandwidths used at the vertices. Any edge whose 
score exceeds a critical value c (c = 0.8 by default) is spht. 
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Table 1. Summary of the parameters used to obtain and fit tfie King models for M62 

Core radius (according to most updated Harris catalog) 0.18 arcmin 

Tidal radius (according to most updated Harris catalog) 8.97 arcmin 

k constant from the King profile (see eq. 11) 12264.7 stars arcmin 

Surface density of non-member field stars from the fit (see eq. 11) 226.4 stars arcmin" 

Surface density of non-member field stars according to Besangon model 216.6 stars arcmin" 
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Table 2. Summary of the parameters used to obtain the Besangon model for the 

non-cluster stars in the M62 region. 



Field of view 

Field 



Extinction law 

Diffuse extinction 
Discrete clouds 

Selection on 

Intervals of magnitude 



Photometric errors 



Small field 
1=353.57°; b=7.32° 
Sohd angle=0.190 square degree 

0.0 mag/kpc 
^^,=1.36; Distance=Opc 

15.75 <B< 25.14 

15.56 <V < 23.71 

14.75 < 7 < 22.44 
Error function: Exponential 
Band=5; A=0.006, B=22.68, C=0.866 
Band=y; A=0.005, B=23.22, C=0.899 
Band=7; A=0.009, B=30.12, C=1.238 



Table 3. Limits for the stars used in our analysis of M62 

V magnitude where the completeness limit is reached 21.51 
AV^ between the completeness limit and the dimmest star in our ground observation 2.2 

Distance where the ratio of GC stars to total number of stars drops to 0.1 6.20 arcmin 
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Fig. 1. — Flow chart describing our technique to map the differential extinction. 
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Fig. 2. — On the left, the density distribution per area of the stars in M 62 as a function of 
distance to the GC center, minus the constant density distribution found for the field (dashed 
green line), is plotted as the solid black line, while the King model for the cluster is plotted 
with a red solid line, and the constant density distribution of stars in the field provided by 
the Besangon model is shown as a blue dashed line (see table [1] for the parameters used to 
build and fit the model). On the right, P{X = l|r), the probability of the stars to belong 
to the cluster as a function of distance from the cluster center. 
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Fig. 3. — On the left, CMD from the Besangon model of the Galactic stars in the field 
centered in the cluster position (see table [2] for the parameters used to build the model). On 
the right, our observed CMD, with the cluster belonging probabilities as a function of color 
and magnitude, P{X = l\c,m), represented by the different colors of the stars as indicated 
in the colorbar. Notice that we restrict our analysis to stars brighter than the completeness 
limit shown in table [31 
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Fig. 4. — Our observed dereddened CMD at four different distances from the cluster center, 
with the probabihty of every star to belong to the cluster as a function of position in the 
sky, color and magnitude, P{X = l\r,c,m), represented by the different colors of the stars 
as indicated in the colorbar. Notice that we restrict our analysis to stars brighter than the 
completeness limit shown in table [3l 
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Fig. 5. — CMD of the example GC, M 62, with the ridgehne provided after the first (upper 
left), second (upper right), and third (lower left) iterations of our method. Warmer (redder) 
colors mean higher density of stars in a region of the CMD. On the lower right, the HB 
ridgeline is also shown, along the different parts of the CMD morphology: main sequence 
(MS), turn-off point (TO), subgiant branch (SGB), red giant branch (RGB) and horizontal 
branch (liB). 



- 39 - 




Fig. 6. — Example of how we move one of M 62 's star along the reddening vector (in green) 
until it intersect the ridgeline. From there a value for the color excess E(B-V) is obtained. 
In the upper left box, we schematically show how the error ellipse moves and how the error 
in the shift is calculated. The 'point of last contact' of the error ellipse with the ridgeline 
is represented in a, the shift of the original datapoint until it intersects the ridgeline is 
represented in b, and the 'point of first contact' of the error ellipse with the ridgeline is 
represented in c. 
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Fig. 7. — Extinction map for the M 62 cluster (left), along with its resolution (middle), and 
its precision (right), as provided by our technique. The x marks the position of the center 
of the cluster. The color code gives the color excesses E{B — V) for the extinction map, the 
bandwidths used in the resolution map, and the standard deviation a of the color excesses 
in the precision maps. Notice that the adopted reddening zero point is defined by where our 
ridgeline lies, and has to be established from other methods (see text). 
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Fig. 8.— B -V vs. V (left) and ^ - / vs. V (right) CMDs of M 62, before and after 
applying our dereddening technique to them. Color bars show the range in the densities of 
stars in the CMD (xlO^ stars per square magnitude). 
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Fig. 9. — Distribution of stars in the B — V color around the calculated ridgeline for two 
different magnitude cuts (1.5 mag above and 0.5 mag below the TO) for the stars in M 62, 
before (blue) and after (red) differentially deredden the CMD. The FWHM of the density 
distributions at the two magnitude cuts narrows by more than a factor of 2 after we have 
differentially dereddened the CMD (for the MS cut, FWHM = 0.262 mag before and 
FWHM = 0.115 mag after; for the RGB cut, FWHM = 0.267 mag before and FWHM = 
0.103 mag after). The profiles also look more unimodal. 
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Fig. 10. — On top, comparison of the extinction map provided 



by our technique (le ft) with 



Schlegel et al. 



fll998l ) (right). 



the extinction map created from the dust temperature maps by 
We have degraded the resolution of our map to make it equal to the Schlegel map. We have 
found an extinction zero point offset between both maps of E{B — V) = 0.47, that we have 
added to our m a p. On the bottom, comparison between the extinction values provided by 



Contreras et al. 



( 120101 ) from the analysis a set of RRLyrae, and the extinction values at the 
same positions provided by our map: on the left, the substraction of both sets of extinction 
values allows to find the extinction zero point for our map, E{B — V) = 0.55, as the weighted 



average of the su 



j stract ed values (green line); on the right, comparisons of both sets of values 



( IContreras et al 



zero point). 



fcOlOl ). as red crosses; this paper, as blue circles, after adding the extinction 



